XUV Opacity of Aluminum between the Cold-Solid to Warm-Plasma Transition 
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We present calculations of the free-free XUV opacity of warm, solid-density aluminum at photon 
energies between the plasma frequency at 15 eV and the L-edge at 73 eV, using both density func- 
tional theory combined with molecular dynamics and a semi-analytical model in the RPA framework 
with the inclusion of local field corrections. As the temperature is increased from room temperature 
to 10 eV, with the ion and electron temperatures equal, we calculate an increase in the opacity in 
the range over which the degree of ionization is constant. The effect is less pronounced if only the 
electron temperature is allowed to increase. The physical significance of these increases is discussed 
in terms of intense XUV-laser matter interactions on both femtosecond and picosecond time-scales. 
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The change in the opacity of material in the XUV 
regime as it makes the transition from a cold solid to 
a warm plasma is a problem that both tests our funda- 
mental understanding of how light interacts with mat- 
ter and is one of increasing experimental interest due 
to the development of ultra-intense XUV free-electron- 
lasers (FELs). For example, the FLASH XUV-FEL 
based in Hamburg now produces 15-fs pulses of XUV 
radiation tunable between 30 and 6 nm, with photon 
numbers > 10^^ per bunch, i.e., peak powers in the GW 
range, at 10 Hz. When focussed with multi-layer-coated 
ofF-axis parabolic mirrors to spot diameters of order a 
few microns, irradiances in excess of 10^^ Wcm~^ can be 
produced - an irradiance regime that has until recently 
remained the province of optical lasers. A thin foil placed 
in the focus of such a beam would be rapidly heated to 
several tens of eV, and during this absorption process 
we expect significant changes in the opacity before the 
onset of any hydrodynamic motion and consequent den- 
sity change. Therefore, to have an understanding of the 
femtosecond heating and provide a predictive capability 
for experimental design, a consistent absorption model is 
required. A complete model, which covers a wide range 
of densities, temperatures, and incident photon energies 
would need to include sophisticated atomic physics that 
accounted for both bound-free and free-free absorption, 
and how these varied with charge state upon heating and 
ionization. The ability to separate the contributions of 
these different mechanisms would provide the opportu- 
nity to test our ability to model the role that each has to 
play in the opacity of dense plasmas. 

It is in this context that the opacity of warm dense Al 
in the XUV regime is of particular interest. With a Fermi 
energy of order 11 eV, plasmon energy just above 15 cV 
and an L-edge at 73 eV, there is a large extent of energy 
space where free-electron theory might be thought to be 
applicable. Indeed, for solid Al, band calculations of the 
density of states (DOS) are in excellent agreement with 
free-electron theory up to the pseudogap between 3d and 
4f bands, at around 40 eV [H . Importantly, calculations 
of the equation of state indicate that for solid density Al 



in thermal equilibrium, ionization of the L-shcU does not 
start to occur until electron temperatures close to the 
Fermi temperature Q. Thus, with a known free-electron 
density and ionic state over a wide temperature range, it 
is an ideal material to test our understanding of free-free 
absorption in warm dense matter (WDM). 

In this rapid communication we present calculations 
for solid density Al at temperatures of up to 10 eV using 
a semi-analytical model for the dielectric response within 
the random phase approximation (RPA) framework and 
with inclusion of particle-hole interactions in terms of a 
local field correction (LFC). We test the validity of the 
model by conducting ah initio molecular dynamics (DFT- 
MD) calculations using the VASP code [| in conjunction 
with a Kubo-Greenwood (KG) calculation of the optical 
properties 0, H. A steady rise in opacity as a function 
of temperature for solid density Al for photon energies 
between the plasma frequency and about 40 eV is cal- 
culated, which is sufHciently large to be experimentally 
verifiable. We conclude by commenting on the remain- 
ing discrepancies between calculated and observed XUV 
opacities for solid-state, room-temperature Al. 

The absorption of light of angular frequency oj by 
a material having a macroscopic dielectric function e^v/ 
is given in terms of the absorption coefficient as k = 
2ci;c^^Im-\/eM [uj). For a quasi- free electron metal the 
macroscopic dielectric function can be expanded in a 
perturbation series of the potential between electrons 
and ions under the assumption that this is small com- 
pared to other energies of the system. Up to second 
order, this expansion has the form (see, e.g., Ref. 0) 
eM{i^) = 1 — Wp/w^ + e'^)((jj). Here ujp denotes the elec- 
tron plasma frequency = nee^/(eome), where rie and 
rUe are the electron density and mass. A systematic ap- 
proach to the calculation of e^^-'(a-'), the first term de- 
scribing the modulation of the free-electron gas due to 
the presence of the ions, was first conducted by Ron and 
Tzoar 71, who's result can be written in the form: 
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Here Vq indicates the Fourier transform of the electron- 
ion interaction potential, S{q) the ion-ion structure factor 
and e{q, w) the electronic microscopic dielectric function. 
This expression has been shown to be exact to second or- 
der in the potential Q . We note that Eq. ([T]) corresponds 
to the Balescu-Lenard approximation in the kinetic the- 
ory approach Further, in the context of the require- 
ment for a theory that is applicable to both the cold-solid 
and hot low-density plasma extremes, we note that it is a 
known result that Eq. ^ is well suited to reproduce the 
classical plasma limit, i.e. inverse bremsstrahlung, when 
the Coulomb potential and classical forms of structure 
factors are considered The applicability of Eq. ^ 
is based on the assumption of quasi free-electron struc- 
ture and can be as such easily extended to a range of 
simple metals. 

Within the RPA, the electron dielectric function can 
be evaluated as a function of temperature throughout 
the degenerate to non-degenerate regime The RPA 
is exact in the high density limit, where long range cor- 
relations dominate. At lower densities, such as those in 
real metals, intermediate and short range correlations be- 
come increasingly important and are described by vertex 
and self-energy corrections. A case of particular interest 
are electron-hole ladder interactions describing excitons, 
which have been shown to give large contributions to the 
absorption strength ^ . We add such contributions to the 
RPA dielectric function in terms of a LFC following the 
work of Higuchi and Yasuhara [l^ . including full tem- 
perature dependence to the RPA polarisation functions 
and calculating the averaged electron screening over a 
finite-temperature Fermi sphere. While this approxima- 
tion might not be sufficient for an in detail study of many- 
body contributions, it does capture the main effect in Al 
which is an increase in the free-free absorption of about 
a factor two in the range of interest. 

In evaluating the ion-ion structure factor in Eq. ([1]), 
there are different approaches that we can use for WDM. 
Numerically, molecular dynamics (MD) techniques 
as well as direct solution of the hyper-netted chain (HNC) 
equations [l^, can be implemented. While accurate 
for ideal classical plasmas, both MD and HNC become 
computationally expensive if the full dynamic response 
of the quantal electron background has to be accounted 
for. In order to simplify the analysis, we use the ap- 
proach initially suggested by Singh and Holz [l^ for a 
liquid metal and applied to the ion-ion structure factor 
of WDM by Gregori et al. [ll] . In this approach, the bare 
ionic response is constructed from the analytical solution 
of the mean spherical equations for a system of charged 
hard spheres embedded in a neutralizing background ■ 
Such an approximation is the equivalent of the Percus- 
Yevick solution for a neutral gas of hard spheres [H, [l^l , 
and has been shown to correctly reproduce MD and HNC 
results at various degrees of inter- ion coupling [l^. The 
idealized one-component plasma (OCP) structure factor 
S'ocp('7) is thus obtained in a fully analytical form. The 
dynamical screening of the electron background on the 
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FIG. 1: (color online) Absorption coefficient versus photon 
energy for solid density Al, Z=3. Comparison between results 
from Eq. ([T]) in the RPA (blue) and with electron-hole LFC 
(purple) and DFT-MD simulations (red). 



ions is then calculated within the linear response frame- 
work SwDuiq) = Socp{q)/[l + Mq)Socp{q)], 

where fv{q) is the attractive screening correction to the 
bare ion-ion interaction [l^, [l^ . The agreement of this 
model with the structure factor computed by DFT-MD 
is seen to be very good. 

The Coulomb interaction potential is not adequate 
to describe the degenerate regime, where the ionic core 
potential is both screened by the electrons and modi- 
fied by the interaction with nearby ions. To take these 
effects into account we use the empty-core pseudopo- 
tential model, where the potential is zero inside a cer- 
tain distance from the ion core, Tc, and —Ze'^/^Aireor) 
outside. Its Fourier transform is given by Vq = 
—Ze"^ cos{qrc)/{eoq^). Here we use Tc = O.6OA, which 
minimizes the difference between the pseudopotential val- 
ues and experimental data obtained through measure- 
ments of the Fermi surface at the reciprocal lattice points 
(111) and (200) of crystalline Al 

We compare the analytic model to DTF-MD calcu- 
lations. In our DFT-MD results, the coupling between 
electrons with different wavevectors is not considered. As 
such, it must be compared to the RPA dielectric model 
without the inclusion of LFC describing electron-hole 
coupling. The DFT-MD absorption was calculated using 
the VASP code along with a KG calculation of the optical 
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FIG. 2: (color online) Absorption coefficient versus tempera- 
ture for solid density Al, Z=3, for 30 eV photons, illustrating 
the absorption peak contribution of the electrons alone (green 
and blue) and of the combined electron-ion system (red). In- 
ner graph shows comparison between results from Eq. ([ij in 
the RPA and with electron-hole LFC and DFT-MD simula- 
tions. 



properties. The ion configurations for the optical proper- 
ties calculations were obtained from snapshots of equili- 
brated molecular dynamics (MD) simulations within the 
framework of finite temperature density functional the- 
ory [2l| . The Al atoms are represented with frozen-core, 
Projector Augmented Wave (PAW) potentials [H, \2^ . 
with the 3s and 3p electrons in the valence band. The ex- 
change and correlation energies are in the Perdew, Burke 
and Ernzerhof form [U of the generalized gradient ap- 
proximation. The DFT-MD calculations were performed 
with the Baldereschi mean value point for sampling the 
Brilliouin zone [1^. For the calculation of the optical 
properties, the Brillouin zone was sampled with the 10 
irreducible fc-points of a Monkhorst-Pack 6x6x6 grid [1^ . 
The plane wave cutoff energy for both the MD runs and 
the KG calculations was 330 eV. The results shown were 
obtained from simulations of 32 Al atoms in a triply pe- 
riodic cubic box corresponding to a density of 2.7 g/cm"^. 
Test calculations with more atoms, higher order /c-point 
sets, and higher plane wave cutoff energies did not signif- 
icantly effect the results for photon energies above 15 eV. 
We have conducted DFT-MD absorption calculations up 
to a photon energy of 45 eV. Absorptions at higher ener- 
gies are not included due to spurious features connected 
to the influence of the pseudopotential on the DOS as the 
L-edge is approached. Gaussian broadening is 0.1 eV. 

The absorption coefhcients predicted by DFT-MD are 
shown in Fig. [1] alongside those obtained within the 
framework given by Ron and Tzoar both in the RPA 
and with LFC. Excellent agreement between the DFT- 
MD and RPA approaches is found in the free electron 
region (up to 35 eV). The discrepancy at higher energies 
is due to the DOS, which is no longer free-elcctron-like, an 



aspect accounted for in DFT-MD but not in the model. 
The LFC result differs significantly from the RPA at low 
temperatures but converges to it for temperatures above 
10 eV, where the electronic LFC contributions become 
small. 

An important aspect of the models is their prediction 
of the opacity as a function of temperature, which we 
show in Fig. [2l for 30 eV photons. Both our analyti- 
cal calculations and the DFT-MD predict a rise in the 
absorption coefficient up to temperatures around 10 eV. 
This absorption rise is due in part to plasmon peak broad- 
ening of the electron liquid and in part to the tempera- 
ture broadening of the peaked, solid-like structure factor 
at low temperatures, as the ions heat up. Both these 
effects enhance the absorption since they allow more k- 
vectors to participate to the process. We note that in 
FEL experiments the ions might not have sufficient time 
to thermalize and only the electron contribution is ex- 
pected. For temperatures higher than the Fermi temper- 
ature, Landau damping is progressively more important 
with a reduction in absorption. 

Whilst some limited information exists about the opac- 
ity of hot dense Al [23] , the error bars on the experimen- 
tal data are sufficiently large to nullify the usefulness of 
comparison with the calculations presented here. Inter- 
estingly, even for solid Al at room temperature significant 
uncertainty in the absorption coefficient still remains. 
There are currently two sets of oft-cited experimental 
data for the cold opacity, given by Henke et al. [1^ and 
GuUikson et al. [29|, which differ by as much as a factor 
two in our range of interest, as can be seen from Fig. [31 
Such differences have large implications in XUV laser and 
high-harmonic generation research, where Al filters are 
often employed. Although more recent measurements us- 
ing XUV lasers conducted at photon energies of 53.7 and 
63.3 eV [lOl show agreement with the Henke data, the 
GuUikson values are still widely used [lH. For the sake 
of completeness, we also present here calculations for the 
cold, sohd case. 

The reduction of the Ron-Tzoar formalism to describe 
absorption in a crystalline solid is straight-forward. In 
this case, the ions form a lattice and the ion-ion struc- 
ture factor S{q) will only depend on those q-vectors 
which match the reciprocal lattice vectors G: <S'(q) = 
<S'soiid(q) — ^TT^m 5(q— G). Inserting this into Eq. ([T|) 
gives the second order term of the macroscopic dielectric 
function for a solid, given in terms of a discrete sum over 
all bands. The same expression was derived by Sturm 
for the study of optical properties of simple metals at 
T = @, ^ . At frequencies above the plasma frequency 
and below the L-shell threshold this is the only significant 
term contributing to the total absorption, the intra-band 
(Drude) term being negligible. 

Figure [3] shows the experimentally obtained absorption 
coefficient for crystalline Al at 300 K of Henke et al., Gul- 
likson et al. and Keenan et al. in comparison with our 
calculations from Eq. ([1]) with S'soiid(q)j with and with- 
out LFCs and with DFT-MD. We note that the RPA 
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FIG. 3: (color online) Absorption coefficient for solid Al. Our 
results with and without electron-hole coupling compared to 
experimental data from Henke et al. [2^, GuUikson et al. [29l ] 
and Keenan et al. [3011. 



is inadequate to reproduce the experimental results, an 
observation in agreement with Sturm, who's calculations 
coincide with our RPA result 



32l |. Good agreement is 



also seen between the RPA calculation and the DFT-MD 
at least up to the pseudo band gap around 40 eV, a result 
that validates the Ron-Tzoar model in the region where 
free-electron theory is applicable. It is immediately clear 
from the plot that the addition of LFC to e(g, lo) presents 
a significant improvement to the calculation, placing it in 
rough agreement with experimental data. This is consis- 
tent with the conclusions of ref. The large discrep- 



ancy between the two data sets however does not allow 
for a more detailed evaluation. That such discrepancies 
exist even in the simple cold solid case suggests the need 
for further accurate experimental investigations. 

In conclusion, we have shown that the consideration of 
appropriate structure factors and potentials can be used 
within the framework given by Ron and Tzoar to model 
the free-free XUV absorption in the cold solid and that 
the extrapolation to the warm plasma state is consistent 
with DFT-MD calculations over a large range of temper- 
atures and photon energies. Furthermore, a rise in the 
absorption as a function of temperature is predicted by 
our calculations, dictated by the behavior of both the 
ion and electron structure factors. Thus experiments de- 
signed to observe the temperature-dependent opacity of 
solid-density Al on different time scales afford the poten- 
tial to substantially increase our understanding of warm 
dense matter. 
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